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ABSTRACT 

We use N-body simulations to test the predictions of the redshift distortion in the 
power spectrum given by the halo model in which the clustering of dark matter par- 
ticles is considered as a result both of the clustering of dark halos in space and of the 
distribution of dark matter particles in individual dark halo. The predicted redshift 
distortion depends sensitively on several model parameters in a way different from 
the real-space power spectrum. An accurate model of the redshift distortion can be 
constructed if the following properties of the halo population are modelled accurately: 
the mass function of dark halos, the velocity dispersion among dark halos, and the 
non-linear nature of halo bias on small scales. The model can be readily applied to 
interpreting the clustering properties and velocity dispersion of different populations 
of galaxies once a cluster-weighted bias (or equivalently an halo occupation num- 
ber model) is specified for the galaxies. Some non-trivial bias features observed from 
redshift surveys of optical galaxies and of IRAS galaxies relative to the standard low- 
density cold dark matter model can be easily explained in the cluster weighted bias 
model. The halo model further indicates that a linear bias can be a good approxima- 
tion only on for k < 0.1hMpc~^. 

Key v^rords: Galaxy clustering - galaxies: distances and redshifts - large-scale struc- 
ture of Universe - cosmology: theory - dark matter 



1 INTRODUCTION 

The power spectrum of the galaxy spatial distribution is 
an important statistic for describing inhomogeneities in the 
Universe. The spatial distribution of galaxies observed with 
a redshift survey is distorted by the peculiar motions of 
galaxies, and a statistically isotropic distribution (e.g. the 
power spectrum or the correlation function) in real space be- 
comes anisotropic in redshift space (Geller & Peebles 1973; 
Davis & Peebles 1983; Bean et al. 1983; Kaiser 1987). On 
the other hand, when measuring clustering for high redshift 
objects in redshift space, choosing the wrong cosmological 
model can lead to an additional anisotropy in the redshift 
space distribution (Matsubara & Suto 1996; Ballinger, Pea- 
cock & Heavens 1996). The theory for the redshift distor- 
tions caused by an assumed world model and by the large 
scale linear motions is now well established. The transform 
of the clustering pattern from the true world model to an 
assumed model is just a simple mapping in the coordinates. 
On large scales and for a linear bias, the redshift power spec- 
trum Pf{k,fi) can be derived (Kaiser 1987) 



(1) 



where ^ is the cosine of the angle between the line-of-sight 
and the k vector, /3 = Q.'^'^ /b, f2 the density parameter, 
6 the linear bias, and Pi^{k) the linear power spectrum of 
dark matter in real space. The hope has been to measure 
the dynamical quantity P and the cosmological parameters 
through studying the anisotropy of galaxies or galaxy clus- 
ter on large enough scales. However, it has been shown that 
the virialized motions within rich clusters (the Finger-of- 
God effect) are so prominent that the clustering pattern 
on large scales of wave number k ~ 0.1/iMpc~^ , is signifi- 
cantly affected (Cole, Fisher & Weinberg 1994; Suto et al. 
1999). The effect of non- linear motions on the redshift dis- 
tortion must be properly modelled in order to measure the 
dynamical quantity and the cosmological parameters from 
the redshift distortion of extragalactic objects. Recent re- 
sults (Peacock et al. 2001) from the redshift distortion of 
the 2dF galaxies further support this point, as even with 
the largest redshift survey available today there exists a tight 
degeneracy in the determinations of the parameters /3 and 
the pairwise velocity (reflecting the non- linear motion). 
As pointed out recently by Jing & Borner (2001, hereafter 
JB2001), the existing analytical model for the nonlinear ve- 
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locity distortion, e.g. the exponential distribution of the rel- 
ative velocity in coordinate space (Davis & Peebles 1983) 
or the Lorentz damping function in k-space (Cole, Fisher & 
Weinberg 1995; Peacock & Dodds 1994), is at best an ap- 
proximation for some scales. Although JB2001 have made 
an extensive study of the redshift distortion for the dark 
matter and for two specific biased tracers in three typical 
cosmological models based on high-resolution simulations, 
it is unknown how to generalize their results to cosmological 
models and/or to biased tracers different from those stud- 
ied in their work without running new simulations. In this 
paper, we present an analytical model for the redshift distor- 
tion based on the halo model, and wc will test the accuracy 
of the analytic model with the results of ,]B2001. 

A key concept in the standard hierarchical scenario of 
structure formation is the formation of dark matter ha- 
los, which are virialized systems of dark matter particles 
formed through non-linear gravitational collapse in the cos- 
mic density field. Since the formation of dark halos involves 
only gravitational physics, accurate analytic models are now 
available for many properties of the halo population, such 
as the mass function (Press & Schechter 1974; Lee & Shan- 
darin 1998; Sheth & Tormen 1999; Sheth, Mo & Tormen 
2001), clustering properties, (Mo & White 1996; Mo, Jing & 
White 1997; Jing 1998, 1999; Sheth & Tormen 1999; Sheth, 
Mo & Tormen 2001; Hamana et al. 2001), and density pro- 
files (Navarro, Prenk & White 1996, 1997; Moore et al. 1999; 
Jing & Suto 2000; Klypin et al. 2001). Such models are very 
useful in understanding the clustering properties of matter 
in the universe, as well as in understanding the bias of the 
galaxy distribution relative to the underlying density field. 
Indeed, since in the hierarchical cosmogony all masses in the 
universe are partitioned in dark halos, the halo clustering 
properties , the halo mass function and the halo density pro- 
files are sufficient for the construction of a clustering model 
for the dark matter in the universe (Scherrer & Bcrtschinger 
1991; Sheth & Jain 1997; Ma & Fry, 2000a,b; Seljak 2000; 
Cooray, Hu & MiraldarEscude 2000; Cooray 2001). Further- 
more, since galaxies are assumed to form through gas cooling 
and condensation in dark halos, a halo model of galaxy clus- 
tering can also be constructed by combining the clustering 
properties of the halo population with an assumption of how 
galaxies populate dark halos (Jing, Mo & Borner 1998; Pea- 
cock & Smith 2000; Seljak 2000; Scocimarro & Sheth 2001; 
Scoccimarro et al. 2001; Berlind & Weinberg 2001). 

The halo model also provides a useful way to under- 
stand the redshift distortion in galaxy clustering (Seljak 
2001; White 2001). In this model, the enhancement of the 
redshift-space power spectrum (relative to that in real space) 
on large scales is assumed to arise in the halo-halo corre- 
lation, while the smearing (Finger-of-God) effect on small 
scales is attributed to the velocity dispersions in dark ha- 
los. In this paper we show, however, that the redshift dis- 
tortion depends on several important effects which were not 
included in early modelling. Using results derived from high- 
resolution N-body simulations wc show that an accurate 
model for the redshift distortion dependents not only on 
the mass function of dark halos, the dark halo clustering 
and the velocity dispersion among dark halos, but also on 
the nonlinear motion of halos at intermediate scales. 



2 POWER SPECTRUM IN REAL SPACE 

In the halo model, the mass density field in the Universe is a 
superposition of halos distributed in space. The mass func- 
tion 4>{M) dM of halos (which is the mean density of halos 
with mass in the range of [M,M + dM]) and their spatial 
distributions n{M, r) dM (which is the number density of 
halos with mass in the range of [M,M + dM] at position r) 
can be described by the extended Press-Schechter formalism 
(Press & Schechter 1974; Bond et al. 1991, Bower 1991, Mo 
& White 1996) or its empirical fitting formula from N-body 
simulations (Jing 1998, 1999; Sheth & Torman 1999). As- 
suming the density profile of halos to be pa (M, r) for mass 
M, the density in the Universe can be expressed as 



p(r)= / n(M, ri)p„(M,|r-ri|)rfAfffri, 



(2) 



Following Peebles (1980) and Scherrer & Bcrtschinger 
(1991), we can divide the 4-dimensional space of r and M 
into infinitesimal volume elements so that the occupation 
number m of halos within any of the volume elements is 
either or 1. Then the density field can be written as. 



p(r) = '^mpaiMi, |r - Til) . 



Using the following relations 



(3) 



{m) = {ni)=<l>{Mi)dMid'vi, 
{ninj)i^j = 4>{Mi)4,{Mj)[l + ^{rij)]d^Tid^rjdMidMj, (4) 

where rij = Ir,— rj |, we can derive the two-point correlation 
function ^(r) for the density field which is a sum of two 

terms. 



(5) 



where the first term on the rhs is the so-called one-halo term 

(.ihir) which is contributed by the internal halo structure 
(the density profile) and the second term is the two-halo 
term £,2h (r) which is contributed by the halo-halo clustering 
(Ma & Fry 2000b): 



P^^ih{r) 



P ^2h{r) 



dMct>{M) p4M, r') pc{M, |r' + r|) , 



■3 / i3 " 

d V d V 



dMi<j}{Mi)pa,{Mi,r') 



j dM2ct>{M2)p4M2,r") 
x^hhiW -r" + rlMi,M2). 



(6) 



In the above equation, p is the mean mass density and 
^hhir, Ml, M2) is the two-point correlation function of halos 
of mass Ml and meiss M2. In k-space, using pa{M, k) to de- 
note the Fourier transform of pa{M,r), where pa(M,k) = 



d^rp„{M,r)e- 
trum P{k), 



we can derive the mass power spec- 



P{k) = Pih{k) + P2h{k) , 

J dM4>{M) pl{M,k), 

jdMi4>{Mi)pc {Mi,k) j dM24>{M2) PciMi 
xPhh{k,Mi,M2) 



P^Pihik) 



P^P2h{k) 
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Where Phh{k, Mi, M2) is the Fourier transform of 
(,hh{r, Ml, M2), that is the cross power spectrum of the halos 
of mass Ml and M2. In the halo model, the mass power spec- 
trum or the mass correlation function is fully determined by 
the halo mass function 0(M), the halo mass density profile 
pa (M, r) and the halo-halo power spectrum. 

The halo mass function 0(M) is given by the Press- 
Schechter formula, 



4>{M) 



•2- 


1/2 


d\na{M) 


.77. 




dlnM 



2/2 



M 



ct(M) 



(8) 



where 5c is a parameter characterizing the linear overdensity 

at the onset of gravitational collapse, which is 1.686 in the 
case of a spherical collapse model in an Einstein-de Sitter 
universe and depends weakly on cosmological models. The 
rms linear mass fluctuation a{M) in spheres of radius R is 
related to the linear power spectrum Pii„{k) 

'ink'^dk 



-f 

Jo 



(277)3 



■PHn{k)W'{kR) 



(9) 



where VK(a::) = 3(sin x — x cos x) /x^ is the Fourier transform 
of a real-space top-hat window function. 
We take the following form 



\i-p 



(10) 



p [r/RsYil + r/Rs 

for the density profile of halos in CDM cosmological models. 
In the above equation, J is a dimensionless density ampli- 
tude, and Rs is a scale radius. The seminal work of Navarro, 
Frenk, & White (1996,1997; NFW) showed that the halos 
formed in CDM models can all be well described by equa- 
tion (10) with p = 1. Subsequent work with higher numerical 
resolution has shown that the innermost density profiles of 
the halos are steeper than the NFW form especially for the 
halos of galactic mass (Fukushige & Makino 1997; Moore et 
al. 1998; Jing & Suto 2000; Fukushige & Makino 2001), and 
the profiles are better described by p = 1.5. Defining the 
radius of a halo, i?2oo, such that within it the mean matter 
density is 200 times the mean density, the density profile is 
characterized by the two parameters p and the concentration 
parameter c = i?2()0 / Rs ■ The concentration parameter was 
found to depend on the halo mass, and can be computed for 
p = 1 semi-analytically (NFW). In this paper, we will use 
the following fitting formula 



c(M) 



c(M) 



7-0 

V M 



P -■ 



3.5 



M, 

"m 



1, 



1.5, 



(11) 



for the LCDM model which fits well the concentration pa- 
rameters measured by Jing & Suto (2000) for this model 
(see also Lokas 2001). A simple fit to the concentration pa- 
rameters of SCDM halos in NFW gives 



c(M) = 10 



M, 

"m 



p = 1.0. 



(12) 



and as for the LCDM, the concentration for p = 1.5 is 0.5 
times of that for p = 1. The mass of a halo M is related 
to i?200 simply by M = ^^RlooP, and M« is defined by 
requiring ^(M*) = 5c- 

The halo distribution is assumed to be linearly biased 
relative to the underlying mass distribution, that is the 




k(hMpc-") 

Figure 1. The real-space power spectrum of dark matter in the 
LCDM and SCDM models. The dotted line is the prediction based 
on the formula of Peacock & Dodds (1996), the solid (for p = 1.0) 
and the short-dashed (for p = 1.5) lines are the predictions of the 
halo model with the density profile cut at the radius R200- The 
long-dashed line is the prediction of the halo model for the density 
profile of p=1.0 without a cutoff in radius 

power spectrum P^hik, Mi, M2) is related to the mass power 
spectrum Pmass{k) 

P^h{k, Ml, M2) = 6(Mi)6(M2) P-mass (k) (13) 

where 6(M) is the bias factor of halos for mass M. We will 
use the analytical formula derived by Mo & White (1996) 



b{M) = l + {i.^- l)/5c 



V = 5c/a(M) , 



(14) 



for the bias factor. Following Ma & Fry (2000b), we will 
use the linear density power spectrum Piin{k) for Pmass{k), 
partly for simplicity but more importantly because it can 
partly account for the exclusion effect between halos. In §4, 
we test the accuracy of this assumption. 

Now, consider two cosmological models. The LCDM 
model is a spatially fiat model with matter density Qm = 0.3 
and cosmological constant Qa = 0.7. Its linear power spec- 
trum has the shape parameter F = Qmh = 0.2 and the am- 
plitude as — 1 (which is the rms mass density fluctuation in 
a top- hat window with radius Sh^^Mpc). The SCDM model 
is a spatially flat model but without a cosmological constant. 
Its linear power spectrum has F = 0.5 and as = 0.6. 

Figure 1 shows a comparison between the real-space 
power spectrum derived from the halo model and the fit- 
ting formula of Peacock & Dodds (1996; PD96) for the two 
cosmological models at redshift z — 0. First, we used the 
fitting formula of Ma & Fry (2000b) for the Fourier trans- 
form of the density profile pc,{M, k). When obtaining the fit- 
ting formula they apparently assumed that the density pro- 
file (equation 10) extends to infinity. With their fitting for- 
mula, we obtained the mass power spectrum for the LCDM 
model ( the long-dashed line in the upper panel) which 
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Figure 2. The Fourier traiisforiii of the density profile pa{c,r). 
The dashed and the dotted Hnes arc for c=l and c=5 respectively 
with the density profile cutoff at radius R200 ■ For comparison, the 
solid line is the fitting formula of Ma & Pry (2001b) in which the 
density profile extends to the infinity 



is in good agreement with the fitting formula of PD96 on 
small scales {k > 2/i~^Mpc), but has much more power 
on large scales. The reason is that the halo density pro- 
files of equation (10) arc assumed to extend to infinity, 
while there should be a radial cutoff in the halo struc- 
ture (Sheth, Hui, Diaferio & Scooccimarro 2001). When 
we take R200 as the radial cutoff for the Fourier transform 
pa{M,k) = jlj^™" d^rpa(M, r)e"* '', we get a mass power 
spectrum which is in very good agreement with the fitting 
formula of PD96 both on small and on large scales in the 
LCDM model. The significance of the cutoff is illustrated in 
Figure 2 for p = 1, where we plot our computed pa{M,k) 
against the algebraic expression without a cutoff of Ma & 
Fry (2000b) (the solid line). The radial cutoff brings about 
a flat Pa{M, k) at small k (large scale), significantly smaller 
than without a radial cutoff for k < l/R2oo- As we have 
seen from Figure 1, the radial cutoff is important to get the 
correct real-space power spectrum, and we will adopt this 
cutoff throughout the paper (A cutoff was also imposed in 
the work of Ma & Fry 2001a,b; CP. Ma, private communi- 
cations). For the SCDM model, the halo model prediction is 
50 percent higher than the fitting formula of PD96 on small 
scales (for k > l/i~^Mpc). A similar discrepancy is found 
for the LCDM model (the figure not plotted in the paper) if 
we replace p (equation 10) with pcrit and take R200 as the 
radius of a sphere of mean interior density 200pcrit, where 
pcrit is the critical density of the universe. Therefore, the 
inaccuracy of the halo model prediction for the mass power 
spectrum is perhaps 50 to 100 percent (in contrast to many 
recent studies which claimed that the halo model is much 
more accurate). In a subsequent paper we will examine this 
issue in a greater detail by comparing the halo model with 
the results from high-resolution N-body simulations of 512^ 
particles. 




1 2 3 4 5 



Figure 3. The 1-D velocity distribution of particles within viri- 
alized halos in simulations (symbols) is well described by a 
Maxwellian distribution with cr^^ /V200 = 1/1.53 (the solid line) 

The mass power spectrum calculated in the halo model 
is quite insensitive to the choice of p for the innermost 
density profile. The difference in P{k) between p = 1 and 
p = 1.5 is rather small for k < 10/iMpc~^, and we will con- 
sider only p = 1.5 for the calculation of the redshift-space 
power spectrum below. We also found that the real-space 
power spectrum in the halo model is quite robust against a 
reasonable change of the mass function 4>{M) and the linear 
bias factor b{M). In contrast we will see that the red-shift 
power spectrum is quite sensitive to the changes of these 
functions. 



3 THE POWER SPECTRUM IN RED SHIFT 
SPACE 

The dark matter halo is nearly virialized inside Ji2oo- The 
velocity distribution of the dark matter within the halo 
should be approximately Maxwellian distributed with a one- 
dimensional velocity dispersion al^ (Sheth 1996). We have 
tested this assumption with the N-body simulations of Jing 
& Suto (1998), and confirmed that the assumption is valid 
(see Figure 3). Furthermore, we found that 

(tI°/V2oo w 1/1.53 (15) 

slightly smaller than l/\/2, where V200 ~ \/ GM/R2m is 
the circular velocity at the viral radius R200, and M is the 
mass of the halo within i?2oo- The virial motion in the halo 
elongates the density distribution of the halo along the line- 
of-sight in redshift space . The density distribution of a halo 
in redshift space pa{T-K, To-) is a convolution of the real-space 
density profile with the velocity distribution along the line- 
of-sight, 

pf(7',r,ra)= / pacir-^ - Vz/H,ra) f(Vz)dVz , (16) 
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Figure 4. The redshift power spectrum of darlc matter (dot lines 
in the upper panels) and contributions from the one-halo and the 
two- halo terms (the lower panels). The dashed lines in the upper 
two panels are the results after accounting for the finite size of the 
simulation box. The crosses arc the simulation results of JB2001. 
In the upper panel, from top to bottom, the values of k are 1.1, 
1.7, 2.7, and 3.1 liMpc^-' respectively, and the values P"{k,fj,) 
are multiplied by 1, lO^"', lO^'^, and lO^"* respectively for clarity. 
For the halo model, the Press-Schechter formula, the formula of 
Mo & White (1996), and Kaiser's formula are used for the mass 
function of halos, the halo bias and the redshift distortion of halos. 
In the lower panel, the one-halo and two-halo contributions are 
drawn by dotted and dashed lines. Their sums are the solid lines. 
The wavenumber is 1.1 and 3.4 hMpc~^ for the upper and lower 
lines and the P^{k, n) values have been multiplied by 1 and 10""^ 
respectively for clarity. 



where r-n and r<T are the distances parallel and perpendicu- 
lar to the Une-of-sight to the halo center, H is the Hubble 

constant, and f{vz) is the velocity distribution along the 
line-of-sight. The redshift space power spectrum P^{k,u) 
can then be written as, 

P^{k,n) = Pih{Kfj.) + P2h{k,n) 

o f T 1.2 2 1D2(m) 

Pfhik,!^) = dMc^{M)pl{M,k)e-'"' 

PUk,l^) = y"dMi<^(Mi)p„(Mi,fe)e-'='^'^"'''(^i)/^ 

X j dM2(t>{M2) Pa{M2,k)e 
X P^h{k,ix,Mi,M2). 



fc^M^o-J" (M2)/2 



(17) 

Where Phh{k, Mi, M2) is the redshift power spectrum for 
two dark matter halos with mass Mi and M2. Under the 
assumption that every halo is moving as a whole according 
to linear theory (Kaiser 1987) , we can derive the redshift 
power spectrum for halos, 

Pihih /i, Ml, M2) = [1 + n° VV6(Mi)] P^h{k, Ml, M2) 

X [1 + n° VV&(M2)] . (18) 



Figure 5. The same as the upper panels of Figure 4 but with the 
modified formulae for the halo mass function and the linear bias 
factor 



We will show that the formula is a good description on large 
scales, but on small scales it gives too much power at large 
p,. A modified formula will be given in the next section. 

In Figure 4 we plot the redshift power spectrum for 
dark matter based on this halo model, and compare it to 
the simulation results of JB2001 for the two cosmological 
models. From the dot line in the top panels, we see that 
the difference in the redshift power spectrum between the 
halo model and the simulations is small on large scales (i.e 
small k). But on small scales the halo model gives signif- 
icantly more power at larger /i than the simulation data, 
especially in the SCDM model. It is also noted that the fi- 
nite box size effect should be taken into account. In Figure 
4 we also plot P^(fc, /x)( the dashed lines in the upper two 
panels) calculated using a cut-off of the power spectrum for 
a scale smaller than the box-size (100/i~'^Mpc in our simu- 
lation). The value of P^(k,ij) at larger /x are lower in the 
case, but the changes axe quite small. The difference is so 
significant that we have to check all the assumptions used 
in our model: the halo mass function; the density profile of 
halos; and the halo-halo redshift power spectrum. Before do- 
ing that, it would be helpful first to check which term, the 
one-halo term or two halo term, is the main contribution to 
the discrepancy. From the lower panels of Figure 4, it can be 
easily found that the two-halo term dominates the redshift 
power spectrum at large p, even for small scales (fc about a 
few /iMpc~^), in contrast to the real-space power spectrum 
(Ma & Fry 2000b). The reason for the difference is that the 
strong real-space clustering power from the one-halo contri- 
bution is suppressed by its internal motion, while the infall 
(merging) of halos to each other enhances the P^{k, p). We 
also find that the small halos contribute significantly to the 
two-halo term. As for the real-space power spectrum, we 
find that the density profile of dark matter has little effect 
on P^{k,n) for k < 4hMpc-^ (the upper limits of k in our 
paper) for 1 < p < 1.5. In the following section, we will check 
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Figure 6. Upper panels - The real-space power spectrum of halos 
with mass M = 1.75 X IO^/i-IMq for the LCDM and M = 
6.0 X W^^H-^Mq for the SCDM model measured from the N- 
body simulations (the symbols) compared with the formula of 
Eq.(13) (the solid lines). Lower panels - the redshift space power 
spectrum of halos in the simulations. The values of k are 1.4, 
1.7 and 2.1 hMpc~^ respectively. The dotted , the dashed and 
the long-dashed lines are the predictions of equation (20) with 
<T = , 200, and 300 km s"'^ and the nonlinear real-space power 
spectrum of halos. The P^^(fc, ^, M) values have been multiplied 
by 1, 10~^,10~'^ respectively for clarity 



the validity of the assumptions about the mass function and 
the halo-halo correlation. 



4 MODIFICATIONS OF THE ASSUMPTIONS 
ABOUT THE HALO DISTRIBUTIONS 

It has been shown by Jing (1998) and Lee & Shandarin 
(1998) that the Press-Schechter formula, equation (8) and 
equation (14) for the halo bias are not good descriptions 
for small halos in high resolution simulations. Actually the 
simulations produce fewer halos than the P-S formula in 
an intermediate mass range 0.2 < M/M* < 1, but more 
halos both for smaller and larger ma sses. Motivated by these 
findings, (sheth & Tormen (1999)) have found the fitting 
formula for the mass function and the halo bias parameter 
(see also Sheth, Mo & Tormen 2001), 



0(M) = 
h{M) = 



p d log v' 
AP dlogM 

t 1 

av — V 

1 + ; 



A 



+ 



1+^(^^^1)1/2. 



//2 



1 + {av'Y 



(19) 



where a = 0.707 , n = 0.3, A = 0.32 and v' = {5^/o{M)f. 
We now take these formula and recalculate the redshift- 
space power spectrum for dark matter. The results are 
shown in Figure 5, where we can find that at small p the 
agreement between the model and the simulation is im- 
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Figure 7. Upper panels — the same as Figure 5, but with the 
nonlinear motions of halos included in the halo model. The pa- 
rameter a is 200kms~l for the LCDM and 300kms~l for the 
SCDM model. Lower panels - the same as the upper panels but 
for galaxies described by the cluster weighted model 



proved very moderately, but not at larger fi. The results 
can be explained as follows: the modified formulae for the 
halo mass function and the halo bias give a higher num- 
ber density and a higher bias factor for small halos (which 
have mass less than 10^^ Mq) than the P-S and Mo & White 
formulae used in the last section, both of which boost the 
two-halo term of the real-space power spectrum. While the 
hnear redshift compression due to the Kaiser effect (which 
increases the redshift power spectrum) is reduced by a larger 
bias, the net effect of using the accurate fitting formula en- 
hances the redshift power spectrum at large fi and at large 
k, thus spoiling the agreement of the halo model with the 
simulation . 

Considering carefully all the assumptions in the halo 
model which may cause the discrepancy with the simula- 
tion results, we find only one assumption still needs to be 
checked, i.e. the redshift-space power spectrum of the ha- 
los. The halos were assumed to be linearly biased relative 
to the linearly evolved matter distribution, with the velocity 
distortion fully described by Kaiser's formula. During the 
non-linear evolution both assumptions about the redshift- 
space power spectrum of the halos may break down. Here 
we check the assumptions with the help of the N-body simu- 
lations of the box size 100/i"^Mpc of Jing & Suto (1998). A 
comparison of the real space power spectrum of the halos in 
the simulations with equation (13) for halo masses around 
i\f = 3 X IO^^Mq (the halos around this mass are most im- 
portant for the two-halo term) are shown in the upper panel 
of Figure 6, where the fitting formula of Jing (1998) is used 
for the linear bias factors b{M). Overall, equation (13) gives 
a reasonable approximation for the real space power spec- 



© 0000 RAS, MNRAS 000, 000-000 



An analytical model for the non-linear redshift- space power spectrum 7 




0.1 1 
k [h Mpo-i] 



0.01 0.1 1 10 
k [h Mpc-i] 



1 1 1 












LCDM 




SCDM 

" 



0.01 0.1 1 10 0.01 0.1 1 10 10= 

k(hMpc-') k(hMpc-') 




Figure 8. Upper panels — the real-space power spectrum of 
galcixies. The left panel is for the optical galaxies of the Las Cam- 
panas Redshift Survey (the data from Jing & Borner 2001) and 
the right panel is for the IRAS galaxies of the PSCz survey ( the 
data from Hamilton & Tegmark 2000). The dashed lines are the 
predictions of the LCDM model with a = 0.08 in the left panel 
and 0.25 in the right panel. Lower panels — the scale dependent 
bias h{k) predicted by the cluster weighted bias model for the two 
CDM models. From top to bottom, the value of a are 0.08, 0.13, 
0.25, 0.32 



tmm of halos in the simulations of both CDM models, with 
the difference less than a factor of 2. More precisely, for the 
small scales k ~ 2hMpc~^ the simulation results are about 
a factor of 2 higher than the adopted formula in the LCDM, 
but the difference between equation (13) and the simulations 
is much smaller in the SCDM model. From these compar- 
isons, we are sure that the differences in the real-space power 
spectrum between the simulations and the adopted simple 
formula cannot fully account for the discrepancy found for 
the redshift space power spectrum in Figure 4. 

The final possible cause for the discrepancy is that the 
redshift distortion of halos is not precisely described by 
Kaiser's formula. Now we check the power spectrum of ha- 
los in redshift space using the simulation data of ,]ing & 
Suto (1998). Some typical examples (the symbols) are shown 
in the lower panels of Figure 6. For comparison, the lin- 
ear prediction of equation (18) is plotted (the dotted lines), 
where the real-space power spectrum from the simulations 
is used for {k, Mi, M2). The linear prediction is much 
higher than the simulation results for larger values of fi. 
Obviously, the redshift distortion of halos is not well de- 
scribed by Kaiser's formula. The result is not surprising, 
as we have noted that the peculiar velocity of massive ha- 
los deviates typically from the prediction of linear theory by 
~ 30% (Colberg et al. 1998). We parameterize the non- linear 
peculiar velocity effect by a Lorentz damping factor. 



Figure 9. The damping factor D{k, fi, a) for the LCDM model. 
The solid lines are our fitting formula (Equation 23) with ay fixed 
to llOOkms-i. The values of a are inserted in each plot 



[i + n°V/6(M2) 



(20) 



P^k,n,Mi,M2) = 



P^-^{k, Ml, M2) r o0-6 ,2/,/.^ 

l + (fcM^.)V2 [^ + " '^/^^'^^ 



where (Jh is a parameter characterizing the pairwise velocity 
dispersion of halos. The parameter may be a function of 
halo mass Mi, M2 and k, but Shcth & Diafcrio (2001) show 
that ah is insensitive to the halo mass. So here we treat 
it as a constant. This is a valid approximation, as shown 
in Figure 6. The dashed lines in the lower left panel and 
the long-dashed lines in the lower right panel agree quite 
well with the simulation results. The corresponding val- 
ues are 200 km s"'^ for the LCDM model and 300 km s"'^ for 
the SCDM model, qualitatively consistent with the findings 
of Colberg et al. (1998). In the upper panel of Figure 7 we 
present the redshift power spectrum for dark matter using 
the three modifications discussed above: the modified Press- 
Schechter function, the modified halo bias (including the 
non-linear bias from the simulation), and the non- linear pe- 
culiar velocity of halos. The result is in very good agreement 
with the simulation result of P* (fc,/i) at all scales (with a 
difference less than a factor of 2). 

In a recent paper, Padilla & Baugh (2002) studied the 
power spectrum of galaxy clusters in redshift space using 
largo N-body simulations. They found that there is no damp- 
ing of the power spectrum at high wavenumber in redshift 
space. It is difficult to make direct comparison between our 
results and theirs because the halos they used in their anal- 
ysis are only massive ones, which are separated by large 
distances and the nonlinear effects may be insignificant. 



5 THE POWER SPECTRUM OF GALAXIES 

This halo model can be readily extended to the real-space 
and redshift-space power spectra for galaxies if we assume a 
model for the number of galaxies N per halo of mass M (or 
the halo occupation number). It was shown by Jing, Mo & 
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Model 


Ai 


71 




A2 


72 




SCDM 


4.43 


0.84 


5.44x10-3 


8.27 


0.90 


1 


LCDM 


8.41 


1.05 


6.75x10"'' 


9.22 


0.84 


0.1 



Table 1. The fitting values of the coefficients in Equation 23 



Borner (1998) that the LCDM model with a power-law halo 
occupation number model, 



JV/Moc (M/Mo)"" 



(21) 



can accurately explain the observed two-point correlation 
function and the pairwise velocity dispersion of the galax- 
ies in the Las Campanas Rcdshift Survey, where a = 
for M < Mo, and a = 0.08 for M > Mo and Mo is 
6 X 10^^/i^^Mq. Since a is positive for M > Mo, i.e. the 
number of galaxies per unit dark matter mass within a mas- 
sive halo decreases with the halo mass, JMB98 called this 
model the cluster- (under) weighted model. They also noted 
that this parameterization is consistent with the observed 
trend of Carlberg et al. (1996) for the CNOC clusters. In 
upper left panel of Figure 8, we compare our result with the 
real space power spectrum of the Las Campanas Redshift 
Survey (LCRS). The dashed line is our result for LCDM 
model and the value of a is 0.08. There is good agreement 
on all the scales. 

It has been recently shown that this model can also 
successfully describe the clustering of the IRAS galaxies in 
the PSCz catalog when a higher a ~ 0.25 is assumed (Jing, 
Borner, & Suto 2001; Scoccimarro & Sheth 2001). The upper 
right panel of Figure 8 shows a comparison of the halo model 
with the real space power spectrum measured by Hamil- 
ton & Tegmark (2000) for the survey. In the figure a value 
of a = 0.25 is assumed, and the model power spectrum 
is in good agreement with the observed result. Although 
the spatial bias of the IRAS galaxies is scale-dependent 
compared with the predictions of CDM models as Hamil- 
ton & Tegmark (2000) emphasized, such non-trivial proper- 
ties of the spatial bias are actually predicted by the cluster 
weighted halo model (sec also Jing et al. 1998). The scale- 
dependent bias fo(fc) = {F{k,a)/P{k,a = 0))^/^ is plotted 
in the lower panels of Figure 8 for several choices of a in 
the two CDM models. The bias factor has a nontrivial de- 
pendence on scale for k > 0.1/iMpc~^. It could be regarded 
as scale-invariant (or a linear bias) only on larger scales. 
Since many observation, like the dynamical measurement of 
(3 based on measuring the peculiar velocity and spatial dis- 
tribution of galaxies, rely on the assumption of a linear bias, 
it is very important to note that such observations must be 
carried out on sufficiently large scales (also see Seljak 2001) 

We also calculate the redshift space power spectrum 
of galaxies for a varying from 0.03 to 0.32. The analyti- 
cal result (the lower panels of Figure 7) agrees pretty well 
with the simulation result of JB2001 for the same cluster- 
weighted model, where a = 0.08 is taken. Following JB2001, 
we consider the damping factor 




Figure 10. The same as Figure 9, but for SCDM model. 



which describes the non-linear motion effect on the redshift 
power spectrum. We have fitted the scaled damping factor 
with the following form 



D{k, /Li, a) 



1 + i/(a) (fc^ia„)2 + 'qg{a)k^ia^Y 



(23) 



which takes into account the fact that the damping func- 
tion is approximately a scaling function of fe/icr„. f{a) and 
g{a.) have the forms e~'*i"^^and e~'*^"^^ respectively. The 
fitting values for the parameters are given in Table 1 for 
the two cosmological models. The effective pairwise velocity 
dispersion is f^^^{a)av, (where cTv is the velocity disper- 
sion corresponding to a = ,ie /(O) = 1), which describes 
how the velocity dispersion changes with the parameter a. 
Also note that f{a) depends on the cosmological model and 
the power spectrum. In our model for dark matter, the ef- 
fective pairwise velocity dispersion decreases by 40% when 
(Ts changes from 1 to 0.7. The formula can also be used to 
compare the theoretical models with the statistic in large 
redshift surveys of galaxies. 

This CLW model has been shown to be qualitatively 
consistent with the semi-analytical models of galaxy formar 
tion which take into account star formation (Sheth et al. 
2001; White 2001). This type of model has been applied to 
CDM models to predict various statistics for the galaxy clus- 
tering(Jing, Mo & Borner 1998; Peacock & Smith 2000; Sel- 
jak 2000; Scoccimarro et al. 2000; Berlind & Weinberg 2001). 
The good agreement shown above means that the halo model 
can be successfully applied to predict the power spectrum of 
various galaxies once their halo occupation number model 
is specified. Thus insight on the galaxy formation can be 
obtained from measuring the clustering of different popular- 
tions of galaxies. 



D{k, n, a) = 



P«(fc,a)(l-|-/?/x^)2 



(22) 
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6 DISCUSSION AND CONCLUSIONS 

In this paper we present an analytical model for the non- 
linear redshift-space power spectrum of dark matter and of 
galaxies based on a halo prescription. The model has three 
important ingredients: the halo mass function, the mass den- 
sity profile of halos, and the halo-halo redshift space power 
spectrum. The predicted redshift power spectrum is found to 
be insensitive to the details of the halo density profiles: the 
NFW density profile and a steeper inner density profile ( Jing 
& Suto 2000) yield an indistinguishable redshift power spec- 
trum. When we use, as many others have done (for the real 
space correlation function of dark matter and for the radial 
averaged redshift power spectrum) the Press-Schechter for- 
mula for the halo mass function, the Mo & White formula for 
the halo bias, and Kaiser's formula for the redshift distortion 
of the halo spatial distribution, we find that the predicted 
redshift-space power spectrum for dark matter is too high at 
large fi to be consistent with the high-resolution simulation 
results of JB2001. The reason why the halo model works 
well for the real space clustering but not for the redshift- 
space power spectrum is that in the latter case the result is 
dominated by the two-halo term on non-linear scales. 

We have checked carefully which assumptions have 
caused the discrepancy between the halo model prediction 
and the simulations. First the fitting formulae of the halo 
mass function and the halo bias from numerical simulations 
(Sheth & Torman 1999) are used in the replacement of the 
analytical formulae based on the Press-Schechter formalism. 
This deteriorates the agreement between the halo model 
prediction and the simulation results , because there are 
more small halos and the spatial clustering of small halos is 
stronger in the modified formulae. We also found that the 
bias of the halo distribution is nearly linear relative to the 
linear density field, so the non-linearity of the bias could 
not be the main contribution to the discrepancy. Instead we 
found that the non-linear motions of the halos, which were 
neglected in previous studies, are the main cause. Once we 
take this effect into account the redshift power spectrum 
based on the halo model agrees very well with the simulation 
results. Furthermore, the redshift space power spectrum can 
be precisely predicted if the halo occupation number model, 
e.g. the cluster weighted model, is given for the galaxies. 

Our results show for the first time that the two-halo 
term can dominate some statistics in the redshift space even 
at small scales. Therefore a halo model based on the density 
profiles of halos and on the redshift distribution of halos 
predicted by the linear theory may become inaccurate for 
some statistics of redshift clustering. For the time being, an 
accurate model (analytical or fitting) for the redshift power 
spectrum of the halos, which includes the effects of the non- 
linear bias and the non-linear motions, is needed for the 
prediction of the redshift space power spectra of galaxies 
and dark matter. 

In combination with the cluster weighted bias model, 
we show that a non-trivial scale-dependent bias is generally 
expected for galaxies in CDM models. The bias could be 
regarded linear only on the scales at the wavelength larger 
than 60/i~^Mpc (also see Seljak 2001). The bias is different 
for different populations of galaxies as observed. The cur- 
rently favored LCDM model can well explain the observed 
features of the spatial bias reported for optical galaxies and 



for IRAS galaxies recently, if the cluster weighted bias model 
is applied. 
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